Quantum control and the Strocchi map 
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g : Abstract 

r~| ' Identifying the real and imaginary parts of wave functions with 

O-i. coordinates and momenta, quantum evolution may be mapped onto a 

^ I classical Hamiltonian system. In addition to the symplectic form, 

■ quantum mechanics also has a positive-definite real inner product 

^ . which provides a geometrical interpretation of the measurement pro- 

cess. Together they endow the quantum Hilbert space with the struc- 
ture of a Kaller manifold. 
^ . Quantum control is discussed in this setting. Quantum time- 

^ I evolution corresponds to smooth Hamiltonian dynamics and measure- 

ments to jumps in the phase space. This adds additional power to 
quantum control, non unitarily controllable systems becoming con- 
trollable by "measurement plus evolution". 

A picture of quantum evolution as Hamiltonian dynamics in a 
classical-like phase-space is the appropriate setting to carry over tech- 
niques from classical to quantum control. This is illustrated by a 
discussion of optimal control and sliding mode techniques. 
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1 Introduction 



The mathematical structures of classical and quantum mechanics are usu- 
ally regarded as essentially different. However, many years ago Strocchi[l], 
by identifying the real and imaginary parts of the wave function with coor- 
dinates and momenta, has shown that quantum evolution may be mapped 
onto a classical Hamiltonian system. In particular this setting suits nicely a 
geometrical interpretation of quantum mechanics. This formulation of quan- 
tum mechanics, first proposed by Strocchi[l], has since been rediscovered 
and extended by many authors [2] [3] [4] [5] [6]. Some applications of the 
Strocchi map to the evolution of finite-dimensional quantum systems were 
considered in [7] and structures relevant to the relation of classical to quan- 
tum mechanics were studied in [8]. A particularly interesting extension to 
the original ideas of the Strocchi map was the realization that, in addition to 
the symplectic form, characteristic of Hamiltonian evolution, quantum me- 
chanics also has a positive-definite real inner product. Together they endow 
the quantum Hilbert space with the structure of a Kaller manifold [3] [6]. 
We will discuss this more general framework of the Strocchi map. However, 
for simplicity and whenever possible, we will adhere to the original intuitive 
coordinate formulation of Strocchi. 

The aim of this work is to discuss quantum control in the geometrical 
setting provided by the Strocchi map. Quantum time-evolution will corre- 
spond to smooth Hamiltonian dynamics in a classical-like phase-space and 
measurements to jumps in the phase-space. This is very different from the 
situation in classical feedback control, where the measurements needed for 
the feedback action are not supposed to change the state of the system or 
to change it only very little. However, in some cases, rather than being a 
nuisance, the state disturbance introduced by quantum measurement adds 
additional power to quantum feedback control and, in particular, it changes 
the question of controllability. Non unitarily controllable systems may be- 
come controllable by "measurement plus evolution" . 

The paper is organized as follows: In section 2 a review is made of the 
properties of the Strocchi map as well as some extensions needed to de- 
scribe quantum evolution and control in both the von Neumann and positive- 
operator-valued measure approaches. In section 3, a new quantum control 
method is proposed which combines free quantum evolution and quantum 
measurement to reach a desirable quantum state. Finally, in section 4, quan- 
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turn analogs of classical nonlinear and optimal control methods are shortly 
discussed. 



2 Geometrical formulation of quantum me- 
chanics. The Strocchi map 

Consider a basis {\k)} in a separable complex Hilbert space H*. A general 
quantum state j'^) is 

= (1) 

k 

Define ^ 

'^k = ^ {Qk + iPk) (2) 

where {qu^Pk} is a numerable set of real phase-space coordinates. Then the 
scalar product in the complex Hilbert space Ti* 

^V) {<l'kQk+PkPk^+i (q'kPk - p'kQk^ " (^','0)} 

k 

(3) 

decomposes into the sum of a positive real inner product 

G{^',i^)=lYl + p'kPk) (4) 

k 

and a symplectic form 

n(^ij',i?j =^J2 {%Pk - Pklk) (5) 

k 

Considering Ti,* — {H, J) as a real Hilbert space H with a complex struc- 
ture J, the triple (J, G, Q) equips H with the structure of a Kahler space 
because 

G(^^',^)=n(i;',Ji?j (6) 

The Schrodinger equation — H lip) becomes the set of Hamilton's 

equations 

d ^ _ d 



dtik = it^ 
d_ 



i.Pk = -T^-H 



associated to the symplectic form {^p' and the "classical" Hamiltonian 

H = ^ ^ {{quqj + PkPj) ReHkj + {pkQj - QkPj) Im-fffcj} (8) 
k,j 

with Hkj = {k\H\j). 

One sees that the time evolution of quantum mechanics is equivalent to 
the classical dynamics of a numerable set of coupled oscillators. What is 
unique to quantum mechanics is the special role played by the symmetric 
form G (^ip ,ipy 

Let S be the Hilbert sphere, that is, the space of normalized {\\ip\\ = 1) 
functions in the Hilbert space H. G (^ip , ■0) defines a metric in S. Consider 
now a measurement of an observable A which, for simplicity, we assume to 
have a (possibly degenerate) discrete spectrum. Let a be a (degenerate) 
eigenvalue of A and Pa the projector on the subspace Va of S associated to 
this eigenvalue. After the measurement of A is performed and the value is 
found to be a , the quantum state changes from ijj E S to ipa = jj^^ ^ ^ ^iid 
the probability to find this value is ||Pa'0||^. The metric G {ip',ip) provides 
a nice geometrical interpretation of the measurement process in quantum 
mechanics. 

Given ip E S and e C 5 it is easy to see that 

(V'-0,V-0) (9) 

is minimal when (p = ipa- Because {ip — — (p) = G {ip — (p,ip — (p) one 
concludes that the measurement projects ip on the element of T4 that is 
closest to Ip in the G— metric. The probabihty for this projection is 

Therefore, whereas the symplectic form Q determines time-evolution, the 
G— metric controls the measurement process. It is the special role played 
by the metric that, in this framework, sets apart quantum from classical 
mechanics. 

In the numerable basis {\k)} of finite or infinite cardinality each 
pure quantum state is represented by a point {if, "jf ) in a "phase-space" 
of dimension 2%. Similarly a mixed state will be described by a density 
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on the same space. The density matrix for a mixed state of the form 
Pi't) = Y,nPn\i^n{t)) {iJn{t)\ {Ylin Pn = 1) becomes, using the notations in 
Eq.(2) 



p{t)= J dtdtp it: t) (^'^ (^) + (^)) (^fe' (^) - W it)) \k) {k' 

(11) 

Using the equations of motion (7) and integration by parts one obtains the 
following classical-like equation for the "phase-space" density 

d _>x dp dm dp dm , 

Y.P it, f) = ■ ^ + ^ ■ ^ = - {p, H} 12 

at oq op op aq 

In the Strocchi map framework, the (unobserved) dynamics of quantum 
states is a continuous symplectic evolution in a phase space. On the other 
hand, measurement of a state is represented by jumps in phase space. Be- 
cause the basis {\k)} is arbitrary, we may suppose that this is a basis of 
eigenstates of the set /C of observables that is being measured. Let the 
state before the measurement be ("^,"p^). When a measurement is per- 
formed and the results registered to be /c, the state jumps from ("g*, p*) to 



if = . efc, = —j£^=eu where el and Cu are unit vectors along the 

A;— coordinate and the A;— momentum. 

In feedback control, the results of a measurement are used to correct 
the driving forces and different corrections will be associated to different 
results of the measurement. Therefore the kind of measurement of interest 
in quantum control is a selective one, that is, one in which the result of 
the measurement is registered. However, in quantum mechanics, even if 
the results of the measurement arc no registered, the state of the system 
(or our information about it) is changed anyway. For such non-selective 
measurements one obtains a mixed state. If, for example the initial state is 
a normalized pure state, corresponding to the phase space vector ("g^, p*'), 
after the measurement the state corresponds to a phase-space density 



^ + p^ / \ •\/^^~+^ 



(13) 
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The above considerations refer to complete quantum mechanical projec- 
tions, that is to (selective or non-selective) outputs of quantum mechanical 
measurements. A description of the behavior of a quantum system under con- 
tinuous observation also exists. It uses generalized quantum measurements 
implemented as positive operator- valued measures. Caves and Milburn[9] 
assume that any measurement takes a certain amount of time and that, in 
the infinitesimal time interval dt, the measurement of the observable A cor- 
responds to the operation 

p ^ PapPa (14) 

with 

(2^)"' 

and f daP]^ (a) Pa (a) = 1. Applying Pa (<^) to any state one obtains a 
superposition of eigenstates of A with eigenvectors centered around a and s 
defines the resolution (or the strength) of the measurement. Notice however, 
that Pa (a) is not a projection. It is supposed[9] to represent a generalized 
selective measurement. For a non-selective measurement, that is, one where 
the results are not recorded, 

p — > y daPA (a) pPa (a) (16) 

Prom Eqs.(15)-(16) and the unconditional evolution ^ — —i [H,p\ , the 
following non- selective continuous observation equation for the density ma- 
trix is obtained 

f = -^[^^,P]-|[A,[A,P]] (17) 

This equation is physically appealing, in the sense that, for example, the 
last double commutator term displays a mechanism for the damping of non- 
diagonal terms in the density matrix. Notice however that the Pa (a) op- 
erations are quantum mechanical projections only in the limit of infinite 
strength, or infinite time. 

Because the choice of basis is arbitrary wc may, without loss of generality, 
choose a basis of eigenstates of the measured observable A. Then, one has 
the following equation for the matrix elements of the density matrix 



d 

Jt' 



p,,. = -z {k\ [H, P]k)-- (A. - \>r p,,, (18) 
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leading to the damping of non-diagonal terms, and A^' being the eigen- 
values of A in the states \k) and |A; ). 

Prom the Strocchi map "phase-space" density point of view, Eq.(18) 
means that, if the initial state is a pure state po (7^) = ^ (7^ ~ ^ {'^ ~ '^), 
(in addition to the Hamiltonian evolution) continuous observation splits the 
density into several components which, when t oo, converge to a density 
as in (13). Independently of the conceptual interest of generalized measure- 
ments and continuous observation, the important point to retain is that the 
operation of measurement modifies Hamiltonian evolution. Hence, it might 
play for quantum control a role similar to the one that is played by dissipation 
in non-linear classical control techniques. 

3 Control by measurement plus evolution 

As seen before, in the Strocchi map phase-space, undisturbed time-evolution 
is a smooth Hamiltonian dynamics in phase-space, whereas measurements 
correspond to jumps in the phase-space. This last aspect is very different 
from the situation in classical feedback control, where the measurements 
needed for the feedback action are not supposed to change the state of the 
system or to change it only very little. However in some cases, rather than 
being a nuisance, the state disturbance introduced by quantum measurement 
adds additional power to quantum feedback control. In particular it changes 
the question of controllability. 

For quantum systems with Hamiltonian 

r 

H(t)^Ho + J2^j(t)Hj (19) 

Huang, Tarn and Clark[10] obtained a general result on controllability in- 
volving the Lie algebras generated by the free and control Hamiltonians. A 
bounded quantum system with finite energy has a finite number N of allowed 
states. For this case a necessary and sufficient condition for controllability is 
that the Lie algebra generated by {HqHi, ■ ■ ■ , Hr} be u (N) [12] (or su (N) 
if phases are not important [11]) because this is the smallest group that acts 
transitively on the complex sphere S^~^. 

When the Strocchi map is used to describe quantum evolution of the 
iV— level system, the phase-space has dimension 2N with coordinates {qk,Pk}- 
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The set of transformations that need to be available to have controllability 
is the U {N) subgroup of SO (2iV) corresponding to the real linear canonical 
transformations, that also preserve the form G (Eq.(4)). It contains 

N {N — 1) /2 one-parameter subgroups of the type 



ql — Qi cos 9 + Qj sin 9 

q'j — — Qi sin 9 + qj cos 9 

Pi — Pi cos 9 + pj sin 9 

Pj — —pi sin 9 + Pj cos 9 



(20) 



N{N-1) /2 of type 



Qi — Qi cos 9 — Pj sin 9 

Qj^ = -Pi sin 6* + g^- cos 6' 

Pi — Pi COS 9 + qj sin 6* 

Pj — qi sin 9 + pj cos 9 



and N of type 



(22) 



q^ = qi cos 9 — Pi sin 9 
Pi = qi sin 9 + Pi cos 9 

Suppose now that A = {HQHi,---,Hr}^^ is a proper subalgebra of 
u{N). Then each orbit of the subgroup G{A) C U (N) does not cover 
S^~^. S^~^ becomes a fiber space with the orbits of G{A) as fibers and base 
U (N) /G{A). Then a goal state ipf can only be reached from i(jq if ipQ and 
ipf belong to the same fiber. The system is not controllable purely by the 
action of the unitary evolution J exp {itH^t)) dr but may be controllable by 
the joint action of measurement plus evolution in the following sense: 

Theorem: Given any goal state ■0/ , there is a family of observables M (ipf) 
such that measurement of one of these observables on any x/jq plus unitary 
evolution leads to ipf if G{A) is either O (N) or Sp (|^)- 

Proof: If G{A) = O (N) or Sp (|-/V) we may choose an orthonormal basis 
{(pi} for S^~^ in the orbit G{A)'iljf. Construct an observable M — Ylii'^iP<t>ii 

being the projector on 0,. Measuring this observable on any state ■0o and 
recording the measured value Ojt the state becomes 0jt and then, by unitary 
evolution, ipf may be reached. 

Remarks: 



8 



(i) Because of both the arbitrary nature of the eigenvalues and of the 
orthonormal basis, there is a large family of observables appropriate for this 
type of control. 

(ii) In the result above the state i^f is fixed, but i/jq is arbitrary. If both 
ipQ and ipf are fixed a much simpler set of controlling interactions Hj may be 
sufficient, ipo being fixed, one constructs the M observable by — 1 vectors 
in the iV — 1-dimensional subspace orthogonal to ipQ plus a single vector in the 
orbit G{A)'iljf , non-orthogonal to V'o- Then, G{A) may be a much smaller 
subgroup than the ones in the theorem. 

One sees that, if properly used, the state disturbing effects of quantum 
measurement, rather than being a nuisance, add controlling power over quan- 
tum processes. 

Examples: 

(i) A simple example of non-controllable 3- level system has been discussed 
by Solomon and Schirmer[13]. Let 

H^Ho + u{t) Hi (23) 

with 

Ho^fx\ OOOl; Hi^d\ 1 1 \ (24) 




Hq and Hi generate the algebra of SO (3), therefore the system is not con- 
trollable by unitary evolution. The Strocchi map evolution equations are 



dt\t \-AQ\p' 



(25) 




A being the matrix I u{t) d u{t) d \ . Eq.(25) splits in block form 



d f q + ip \ _ . f -A \ I q + ip ^ ^^^^ 



dt \ q — ip J \ A J y q — ip 



exhibiting the SO (3) nature of the control. In terms of the Hilbert space 
wave functions, Eq.(26) means that ■0* cannot be reached from tp. Three 
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one-parameters subgroups of the SO (3) control group, to be used later on, 
are 

(0) : 






cos 6 


^sin^ 


^sin^ 



I (cos 9 




1) 



;cos^ + 1) 


^sin^ 




-isin^ 



i(cos^ + l) 



i(cos^- 1) 



-isin^ 



-isin^ 


cos^ 






^sin^ 


(cos 6' 


(cos^ + 1) y 



1) 



h2 (e) : 




cos9 + 1) 



sm I 



cos 6*) 




-isin^ 

cos^ 
^sin^ 









i(l -cos^) 



^sin^^ 
(cos^^ 









1) 



1 







COS 0+1) 
sm.9 
- COS 9) 







^sin^ 

cos^ 
^sin^ 



/ cos^ sine \ 

1 

cos^ -sine 

-sine cose 

1 

\ sine cose J 



\ 







(1-cose) 

^sine 

(cose + 1) / 



(29) 



Although not controllable by unitary evolution alone, the system is con- 
trollable by measurement plus evolution. Consider, for example. 
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ipf = j-g" = (0, 0, 0) , "p" = (^^, 0. ^ j I as the goal state. By applying 
hi (— l) and /i2 (f ) hi (— f ) to this state one obtains an orthogonal set 

V'i = {Y = (0, 1,0), y = (0,0,0)} 

V'2 = |^=(^,0,^),^= (0,0,0)1 (30) 

ip3 = [t = (0,0,0), t = (^,0,^)1 

Denoting by Pi the projectors on the states ijji, measurement of an arbitrary 
state by any one of the observables in the family 

^a,P, (31) 

i 

projects it on the fiber of -0/ and then, by unitary evolution, ipf may be 
reached. 

(ii) So far, control by measurement plus evolution has been discussed 
for finite-dimensional spaces. However the same technique may be used in 
infinite-dimensional spaces to reach a large number of states. This is il- 
lustrated for kicked motions in the torus. Let 1^ = {xi,X2^ [— 7r,7r)) be 
coordinates in the 2-torus and the system Hamiltonian H be 

H = - Y+Xl I ^1*1 (t) ( • A • iV iV • >1 • x") (t) xi + U3 (t) | 5 (i - nr) 

n ^ ^ 

(32) 

The switching functions Ui it) take values or ±1 and the matrix A is chosen 
such that M = exp (A) is a hyperbolic 2x2 matrix with integers entries and 
determinant one, this being the condition that insures unitarity of the FIo- 
quet operator [16]. The system is a controlled version of the configurational 
quantum cat [14], a system that describes a charged particle acted upon by 
electromagnetic pulses. When ui (t) ^ the Floquet operator has continu- 
ous spectrum and quantum chaos, in the sense of positive quantum Lyapunov 
exponents[15] [16]. The free and kicked components of the Floquet operator 
are 

Uo = exp fif r) 

Ui = explicit • A»iV + iV • A»^)) , . 

U2 — exp {—ixi) 
U3 = exp {-ix2) 
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Uo corresponds to free propagation, Ui to the action of a linear vector field 
and U2, U-i to scalar potentials. The eigenstates of momentum form a numer- 
able normalized basis 



which is dense on the Hilbert space of the system. The kicks U-i act on these 
states as follows 



By the action of kicks of type U\ the momentum eigenstates move along 
hyperbolas. The kicks U2 and move between different hyperbola. In 
between kicks, free propagation just changes the phase of the states. 

Given now an arbitrary state ip , measuring its momentum the state 
becomes projected on a momentum eigenstate with known eigenvalue, if the 
result of the measurement is recorded. By switching on the appropriate 
sequence of kicks Ui it is then possible to reach any momentum eigenstate. 
Therefore one sees that with measurement and three controlling fields one 
can, from an arbitrary ip, reach any state in an infinite-dimensional dense 
set. 

4 Nonlinear and optimal control 

Classical control is a very mature field where many useful techniques and 

results have been found, many of them still without parallel in quantum con- 
trol. The Strocchi map, yielding a picture of quantum evolution as Hamil- 
tonian dynamics in a classical-like phase-space, may be the appropriate tool 
to carry over techniques from classical to quantum control. We will give two 
examples: 

(i) Optimal control 

Optimal control is an important issue both in classical and quantum 
control and, in quantum control, it has been discussed using variational 
techniques[17] [18] [19]. However, in classical control, Pontryagin maximum 
principle[20] provides a more general framework in the sense that it does not 
require differentiability and can handle piecewise continuous and magnitude- 
limited control. We now show how to carry over this principle to quantum 
control using the Strocchi map. 




(34) 




(35) 
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In addition to the 2A'" phase-space variables {xk = qk',Xk+N = Pk) with 
dynamical laws 



_ <i„ _ 9 M - 



obtained from (7), we introduce a variable Xq with dynamical law 

^ xo = Xo (x, u, t) (37) 



with 

l-T 



I Xo{x,u,t)dt (38) 
Jo 



being the performance functional to be minimized. If, for example, minimal 
controlling energy is desired Xq {x, u, t) = \u{t)\ , etc. 

Then, for each variable in the set "2^ = {xq, qi, - ■ ■ , qN,Pi, ■ ■ ■ ,Pn) an ad- 
joint variable 4>i is defined and a new "Hamiltonian" H [x, (p, u) 



2N 



H {x, 4>,u) = ^ (piXi {x, u) (39) 



i=0 



For each specified initial {xu = qk = Q'fc(O), x^+n = Pk = Pfe(0))and final state 
(xfc = qk = qk{T),Xk+N = Pk = Pk{T)), the optimal control that, in time T, 
minimizes the functional F is obtained by integration of 



'^^k = 4-H(a;,0,M) 



dt^k = -^H(x,0,«) (40) 



it — arg max-yi>gj^H (x, 0, u) 



VL being the domain of allowed controls. 

In the variational formulation of optimal quantum control, the technique 
that has been used is to minimize the functional 

J = (V'(T)|(l-n)|V'(T))+ / Xoi^P, u,t)dt+ [ Im{C {t)\idt-H\^{t)) 

Jo Jo 

n being the projector onto the target state. The first two equations in (40) 
correspond to the Schrodinger equation and to the equation for the Lagrange 
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multiplier ( (t) , whereas the third one corresponds to the equation obtained 
by variation 6lt of the control parameters. However the set (40) is more 
general in that it does not require differentiability in of J and allows the 
specification of arbitrary control domains Q. 

(ii) Sliding mode techniques in quantum control 

A very robust tool used in classical control is the technique of variable 
structure control leading to shding modes[21]. The design of a variable struc- 
ture control has two steps. First, a switching surface must be chosen so that 
the dynamical system restricted to the surface has the desired dynamics. 
Second, a switched control must be found to drive the system to the switch- 
ing surface and, upon interception, to maintain it there. For this step a 
Lyapunov function approach is used, the gradient of the Lyapunov function 
being negative in the neighborhood of the switching surface. In this way 
the tangent vectors to the state trajectory point towards the surface. The 
system is attracted to the switching surface and, once having intercepted it, 
remains there for all subsequent times. Then, the state trajectory is said to 
be a sliding mode. 

From the very nature of the technique one sees that, at least close to 
the switching surface, the dynamics must have a dissipative component and 
therefore a purely Hamiltonian control cannot be used. However, as seen in 
Section 1, observation introduces non-Hamiltonian effects which might play 
for quantum dynamics the same role as dissipation in classical control. 

Consider again the non-controllable example of the previous section and 
suppose that one wants to stabilize the middle energy level. Starting from 
an arbitrary initial state one measures the energy. If the result is zero we are 
done because it means that the state was projected on the middle level. Oth- 
erwise if the result is +11 or —/i, one knows from the controlling subgroups 
in Eqs. (27-29), that there is a control that changes the upper or the lower 
level into a state with a 50% probability of being projected on the middle 
energy state by an energy measurement. One performs this control plus mea- 
surement operation until the result of the measurement is zero. Afterwards 
the control is switched off. If there are no disturbances the state remains 
in the middle level. Otherwise if there are some decoherent interaction with 
the environment, the state should be periodically measured. If the distur- 
bance is small or if the intervals between measurements are small, there is a 
high probability that the system will be projected back on the middle level. 
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without any need for further controlhng operations. 



5 Conclusions 

In addition to a formulation of quantum control in a symplectic geometry 
setting, the main result of this work is the proposal of a new protocol for 
quantum control, which we have called "control by measurement plus evolu- 
tion" . It extends the scope of quantum controllability and is applicable both 
to finite and infinite-dimensional level systems. 

Different aspects of quantum control are united in the framework of the 
Strocchi map, allowing new insights in optimal control and nonlinear control 
techniques. We think that the potential of this picture is not exhausted and 
more analogies with classical problems may be used to obtain progress in the 
quantum domain. 
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